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THEORETICAL CALCULATIONS OF SUPERSONIC WAVE DRAG AT 

ZERO LIFT FCR A PARTICULAR STORE ARRANGEMENT 

By Kenneth MargoUs, Frank S. Malvestuto, Jr., 
and Peter J. Maxie, Jr. 

SUMMARY 


An analysis, based on the linearized thin-airfoil theory for super- 
sonic ’speeds, of the wave drag at zero lif t has been carried out for a 
simple two-body arrangement consisting of two wedgelike surfaces, each 
with a rhombic lateral cross section and emanating from, a common apex. 

Such an arrangement could be used as two stores, either embedded within 
or mounted below a wing, or as auxiliary bodies wherein the upper halves 
could be used as stores and the lower halves for bomb or missile purposes. 
The complete range of supersonic Mach numbers has been considered and it 
was found that by orienting the axes of the bodies relative to each other 
a given volume may be redistributed in a manner which enables the wave 
drag to be reduced within the lower supersonic speed range (where the 
leading edge is substantially subsonic). At the higher Mach numbers, the 
wave drag is always increased. If, in addition to a constant volume, a 
given maximum thickness -chord ratio is imposed, then canting the two sur- 
faces results in higher wave drag at all Mach numbers. For purposes of 
comparison, analogous drag calculations for the case of two parallel 
winglike bodies with the same cross-sectional shapes as the canted con- 
figuration have been included. Consideration is also given to the favor- 
able (dragwise) interference pressures acting on the blunt bases of both 
arrangements . 


INTRODUCTION 


The magnitude of the supersonic wave drag of complete airplane con- 
figurations is known to be dependent not only on the direct drag effects 
generated by each component part of the configuration but also on the 
i ndirect or interference effects introduced by each component on all 
other components. For airplanes that require external fuel tanks, promi- 
nent nacelles, or other auxiliary bodies for storage or bomb and missile 
purposes, the location of such bodies relative to each other and to other 
airplane components has been shown to be an important consideration from 
the standpoint of obtaining low drag. In this connection, reference 1 
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points out that judicious positioning of the components of a multibody 
arrangement can give rise to beneficial interference effects that are of 
sufficient magnitude to allow a volume increase of 25 percent and at the 
same tine actually decrease the wave drag at zero lift. Although the 
possible increases in other forms of drag (e.g., friction drag, base drag, 
and drag due to lift) might partially cancel or even completely nullify 
such reduction in zero -lift wave drag, the importance of positioning 
auxiliary bodies appears clear. 

Inasmuch as interest has been evidenced in the lift and side force of 
conical stores oriented relative to each other, it was believed worthwhile 
to consider also the drag problem for a similar arrangement. In order to 
allow a rigorous analytical solution as well as to simplify the calcula- 
tions, the present paper treats the case of two simple wedge -type stores, 
canted relative to each other, each with a rhombic profile in the span- 
wise direction and both emanating from a single apex. This simple two- 
body arrangement could conceivably be used as two stores, either embedded 
within a wing or mounted below a wing, or as auxiliary bodies wherein the 
upper halves could be used as stores and the lower halves for bomb or 
missile purposes. The zero-lift wave drag of the system is evaluated 
herein by means of linearized supersonic flow theory for the entire super- 
sonic speed range. Corresponding calculations for the case of two parallel 
winglike surfaces with the same cross-sectional shapes as the canted 
arrangement axe included. Consideration is also given to the pressures 
acting on the blunt bases of both configurations. 


LIST OF SYMBOLS 


x,y,z Cartesian coordinates 

z-g z -coordinate defining airfoil section 

g,rj Cartesian coordinates of source points 

V free -stream or flight velocity 

p density of air 

Ap pressure increment 

1 P 

q. dynamic pressure, ipV 


* 


* 


M 


Mach number 
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c max 

"^max 

S,A 


Sy^ Ay 



c D l 


R 

V ] 

Z 
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Mach number parameter, y M^ - 1 

disturbance-velocity potential 

slope of surface (taken in stream direction) 

leading-edge sweepback of body (figs. 1 and J) 

leading-edge sweepback of wing in which canted arrangement 
is embedded (fig. l) 

inclination of inner edge for the canted arrangement 
(positive) 

slope of ridge line, ^(tan 8 + cot A) 
maximum chord 


maximum thickness 


plan-form area and aspect ratio, respectively (For the 
canted arrangement, the area includes the cutout inner 
portion but excludes any wing area exterior to the body- 
leading edge so that S = (c m»y )^cot A and A = 4 cot A 

(fig- D) 

plan-form area and aspect ratio, respectively, of a delta 
wing in which the canted arrangement is assumed to be 
embedded (Sy = (c may ) cot Ay and Ay = 4 cot Ay (fig. l)) 


wave-drag coefficient, 
wave-drag coefficient. 


Wave drag 

Wave drag 
<lSw 


region of integration 
volume of body 


distance between apexes for the parallel arrangement 

semispan of one body for the parallel arrangement 

functions used in appendix A for summarizing expressions 
for velocity potentials 


* * ,F 8 
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Subscripts : 

1,2 used as subscripts on A to denote specific slopes 

1,2, 3, 4 used as subscripts on 0 to denote the fpur basic 

potentials for the canted arrangement 

1A,1B,1C,2A,2B used as subscripts on 0 to denote special forms 

of 0^ and 0 2 

base int base interference 


ANALYSIS 


The analysis is based on supersonic thin-airfoil theory and on the 
assumptions of small disturbances and a constant velocity of sound 
throughout the fluid. These assumptions lead to the linearized equation 
for the velocity potential 0: 

(1 - M 2 )^ (!) 


where M is the Mach number of the flow and the derivatives are taken 

with respect to the variables x, y, and z of the Cartesian coordinate 

system. The general expression for the linearized perturbation -velocity 
potential in space due to a distribution of source and sink singularities 
in the z = 0 plane is (see refs. 2 and 3) 

^z) . - x rr _ , = = ( 2) 

uu R ^(x - |) 2 - p 2 (y - Ti) 2 - |3 2 z 2 

where x, y, and z are the coordinates of the field point (that is, 

the point at which the potential is desired) and | and t] are the 

coordinates (analogous to x and y) of the singularities. The func- 
tion A(|,Tj) represents the particular distribution of singularities 
and is, of course, dependent upon the boundary conditions imposed. For 
the case of wing thickness distributions that are amenable to thln- 
alrf oil-theory calculations, the source -sink distribution function is 
related to the particular thickness distribution involved and is given 
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Jz B“° 


(3) 


Hie integration indicated in equation (2) is performed over the region R 
that is enclosed by the traces in the z = 0 plane of the Mach forecone 
emanating from the point (x,y,z) and by the wing plan-form boundaries. 
Inasmuch as the present problem is concerned with the drag and hence the 
surface conditions, consistent linearization also requires the potential 
to be evaluated in the z = 0 plane. ' Thus, the applicable form of 
equation (2) becomes 


0(x,y) = - 1 




)/(* - O 2 - P 2 (y - Tj) 2 


d| dTj 


00 


The parameter {3 equals i/m 2 - 1 and its reciprocal l/p is the slope 
(absolute value) of the Mach lines. 

Calculation of the wave -drag coefficient involves an integration of 
the tangential components of surface pressure which are in turn related 
to the x-derivative of the velocity potential; that is. 


% 


Si 


C D = I II tr- A(x,y)dx dy 


S.l.j s a 


* - S A(x ' y)ta ^ 


(5) 


where S is the plan-form area and the quantity — = - % 

q V dx 


The particular thickness distribution in the present problem allows 
an additional simplification. Inasmuch as the quantity A(x,y) is a 
constant value over any one region of a wedge -type surface, equation (5) 
may be expressed simply in terms of the potential by performing an inte- 
gration with respect to x; that is, 

CD ’-wf b X0(x ’ y> 


T.E. 

<3y (6) 


ft 


L.E. 
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where the symbol b denotes that the y-lntegration is performed over 
the span. It should be noted, however, that the expression for A0(x,y) 
will generally be different for each wedge region and, therefore, the 
evaluation between the leading edge and trailing edge as required in the 
integrand of equation (6) must be taken in finite steps. The solution to 
the problem consists primarily of calculating the potential for the vari- 
ous surface regions by use of equation (4) and then calculating the drag 
coefficient by means of equation (6). 


Consider now the wing -type arrangement in detail (fig. l). Four 
basic potentials and 04 are required in order to deter- 

mine the wave drag for the entire supersonic range of Mach numbers; the 
areas of integration for each of these are indicated in figure 2. In 
addition, and assume different forms for various conditions. A 


complete summary of the cases treated, mathematical conditions involved, 
and the appropriate forms of the potentials is given in table I. Before 
the potentials are evaluated as indicated by equation (4), appropriate 


expressions for the slopes 


|r z B^n) 


must be obtained. For the out- 


board panels, the slope is found to be 






(7) 


and, for the inboard panels. 


3 „ 1 _ (tan A) (tan 8)t max / c raax 

B ^ 1 - (tan A) (tan 9) 


The limits of integration imposed by region R are readily obtained 
from the information given in figure 2. Expressions for each of the 
potentials obtained as a result of performing the mathematical operations 
indicated in equation (4) are presented in appendix A. In accordance 
with equation (6) the wave -drag coefficient may now be calculated; the 
resultant expressions are rather lengthy and are presented in appendix B. 
As a convenience to the reader and as an aid in performing other calcula- 
tions Involving similar integrals, evaluations of some specific indefinite 
integrals that appear frequently in the analysis herein are given in 
appendix C . 


V 
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Wave -drag calculations analogous to those for the canted arrangement 
were carried out for the case of two parallel bodies separated from each 
other by an arbitrary distance l (see fig. 3). T he mathematical details 
of obtaining the potentials and so forth are omitted herein; only the 
required potential expressions and final drag equations are presented in 
appendix D. It might be noted, however, that the procedures used to 
obtain these expressions are identical to those utilized for the apex- 
adjoining arrangement, although much less effort is involved. For 
example, only two basic potentials are required and the slopes of both 
the inboard and outboard panels are the same. Furthermore, the work of 
reference 3 can be used to determine the potentials and wave drag at all 
supersonic Mach numbers for the case where no interference effects are 
present. 


RESULTS AND DISCUSSION 


As indicated in the analysis, closed-form formulas are presented 
in the appendixes for the complete Mach number range. These formulas 
enable the drag coefficient for any given plan-form arrangement to be 
obtained readily by straightforward numerical computation . Therefore, 
rather than attempting to summarize all possible cases by means of 
plotting series of design charts, the present paper emphasizes only those 
salient points that are borne out by the detailed calculations. 

The obvious question is whether there is any possible advantage drag- 
wise in canting the stores or winglike surfaces as indicated. The intro- 
duction of a finite value of 0, that is, opening the two panels of the 
twin -wedge arrangement relative to one another, gives rise to two effects 
which influence the resultant wave drag. From purely geometric considera- 
tions, it is seen that for a given leading -edge sweepback and a constant 
thickness -chord ratio an increase in the angle 0 results in an increase 
of surface slopes in the stream direction. The trend toward higher slopes 
is even more pronounced if, instead, the thickness -chord ratio is allowed 
to vary and the condition of constant volume is imposed. Inasmuch as the 
supersonic wave drag at zero lift is dependent on the square of the magni- 
tude of such slopes, the effect on the drag is adverse; that is, the drag 
is increased at all Mach numbers. Canting the two panels, however, also 
changes the pressure distribution and, since both positive and negative 
pressures are induced in the field, an assessment of this interference 
effect is required to determine the net change in the wave-drag coefficient. 

Consider the case of 9=0 wherein the arrangement becomes simply a 
delta wing with a lateral cross section composed of two diamonds. For this 
limiting case, the importance and magnitude of the increased-slope effect 
may be assessed directly by means of comparing the drag with that obtained 
for the more conventional wing with a lateral cross section composed of a 
single diamond (equations given in appendix B) . The wave -drag coefficients 
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for both cases are plotted against the aspect-ratio — Mach number param- 
eter A(B in figure 4. The two configurations have, for a given value of 
maximum thickness -chord ratio, the same volume and longitudinal distribu- 
tion of cross-sectional area and the same maximum thickness at each cross 
section. The main difference is, of course, in the distribution of thick- 
ness. The comparison indicates two results which are to be expected in 
view of previously published wave -dreg analyses: First, higher drag 

values are obtained for the "twin-wedge" case throughout the supersonic 
speed range (the effective slopes in the stream direction are much higher 
for this case and hence the drag is correspondingly higher) and, second, 
this increased-drag effect is more pronounced at the higher Mach numbers 
(where the leading edge is supersonic). As an interesting sidelight, it 
might be noted that, inasmuch as the longitudinal distributions of cross- 
sectional areas are the same for both configurations, application of the 
"transonic area rule" would yield the same drag coefficient for both 
cases. Thus, the differences in the ordinates of the two curves give a 
direct measurement of the degree of inaccuracy involved when this rule is 
applied to the supersonic speed range. 

As might be suggested from consideration of figure 4, the lower 
supersonic range (where the leading edge is substantially subsonic) 
appears to be the logical range in which to anticipate net drag reduction 
as a result of canting the two surfaces . Detailed calculations based on 
the equations given in appendix B and covering the complete supersonic 
speed range substantiate this expectation - only in the lower supersonic 
speed range and for relatively small opening angles is the interference 
effect sufficiently large to overbalance the adverse drag effect due to 
the higher local streamwise slopes. At the higher Mach numbers where the 
leading edge approaches the sonic condition or is supersonic, the wave 
drag for the canted arrangement is always higher. Figure 5 presents some 
of the calculations which illustrate these results for the lower super- 
sonic speed range and includes for comparison purposes the uncanted 
arrangement (that is, 3 tan 0=0) previously shown in figure 4. 


The results of figure 5 are directly applicable to cases for which 
the maximum thickness -chord ratio is constant. The curves may be 
replotted in terms of constant volume or for any other geometric con- 
sideration by introducing the appropriate multiplicative factors for 
each point. Consider, for example, the condition of constant volume; 
the volume V B enclosed by both surfaces is given by the expression 

V B . W ( C|m *) (cot A - taa 9 ) 



cot 

°max 


A - 


p 


f 


► 


or 


3 tan 9 ) 


(9) 
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and the maximum thickness- chord ratio by 


3PTb/(w) 5 

Cmax Mot A - p tan 6 


max 


( 10 ) 


The wave-drag coefficient given in figure 5 Is based on the plan- 
form. area S (which includes the cutout portion) and is represented as 




' fc max 
\ c max 


2 

f (p cot A, p tan 8 ) 


(H) 


It is readily apparent from equation (9) that, for a fixed overall length, 
the volume may be maintained constant by either adjusting any two or all 

three of the parameters A. 8. and Thus, for a given volume 

°max 

condition, it is desirable to base the wave-drag coefficient on an area 
which is completely independent of the canted-body geometry so that a 
true indication of the drag itself may be obtained. A corresponding 
drag coefficient Cp’ baaed on the area Sy- may then be written as 


^0Aw) PCpi 1 



2 

(Ij-p cot A)pCp 


( 12 ) 


where £v and Aw are the area and aspect ratio, respectively, of a 
delta wing in which the canted body is assumed to be embedded. For given 
values of the parameters p tan 8 and p cot A, equation (10) gives the 
specific thickness ratio required to maintain the desired volume. Equa- 
tion (ll) in conjunction with figure 5 determines the drag coeffi- 
cient pC D , and then equation (12) may be used to plot results for the 

constant-volume condition. 

Figure 6 presents results obtained for a specific value of the 
3f% 

volume parameter -=• = 0.02. As previously discussed, the ordi- 

( c max)^ 

nate (pAwjfpCD’) of figure 6 gives a direct measurement of the drag 
itself since A^ and ^area used for nondimensionalizlng C-q') are 

independent of the canted-body geometry and, therefore, any plotted point 
in figure 6 may be legitimately compared with any other plotted point to 
determine whether the drag is decreased or increased when geometry param- 
eters are changed in such a way as to maintain constant volume. 

The conclusions to be drawn from figure 6 are as follows: For a 
given overall length and a given Mach number (constant p). It Is possible 
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to redistribute a given volume in a manner practical for store purposes 
and obtain drag reduction by canting the bodies slightly. (Compare, 
for example, the values indicated by the filled- in circles.) If, in 
addition to a constant volume a constant maximum thickness-chord ratio 
is desired, then canting the bodies will always result in increased drag 
(follow the dashed lines) . Similarly, if in addition to constant volume 
a constant sweepback angle A is desired, then canting the bodies will 
result in increased drag (follow vertical lines) . Figure 6 is, of 


course, directly applicable to the volume parameter — — = 0 . 02 , but 

( c max)^ 

the results are indicative of those found for other values of the volume 
parameter. Figure 5 in- conjunction with equations (10), (ll), and (12) 
may be used as just outlined to obtain detailed curves for other volume 
conditions and thickness ratios. - 


Calculations of the supersonic wave drag for the parallel -body 
arrangement consisting of two surfaces of delta plan form, each with a 
simple wedge profile and with parallel axes of symmetry (see fig. 3 ) > 
are presented in figure 7 . The interesting point to be noted is that 
when the two surfaces are in the interference fields of each other the 
resulting interference drag is additive; that is, the drag of the 
arrangement is least when the combination of Mach number and lateral 
displacement of the two apexes is such that the disturbance field of one 
surface does not influence the drag of the other surface. This result 
is expected inasmuch as the pressure due to a single wedge 1 b the same ^ 
in sign over the wedge surface and in the field beyond the wedge; there- 
fore, the introduction of a similar wedge into the field (such as in the 
present case) will result in additional drag of the same sign as the * 

pressure drag of the original body. 

Inasmuch as both the canted and parallel arrangements treated In the 
present paper have blunt bases, it is advisable to point out that the drag 
calculations discussed thus far do not take into account either the basic 
base drag or the Interference -drag contribution resulting from the pres- 
sure field generated by one panel acting on the blunt base of the opposite 
panel. Actually, this interference effect gives rise to a negative drag, 
or thrust, and could conceivably be of the same order of magnitude as the 
interference contribution previously considered. 

In order to assess the importance of the base -drag-interference con- 
tribution, the interference pressures acting on the base have been derived 
for both arrangements of bodies and are presented in appendix E. Specifi- 
cally, the formulas give the interference -pressure coefficient (— ) 

' ^ * base int 

acting along the base line In the plane_of symmetry, that is, x = c^ ny and 
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z = 0. The formulas were obtained by utilizing equation (4) to find the 
interference velocity potential (changing the region R to include only 
that portion bounded by the Mach trace and plan-form boundary of the oppo- 
site or interfering body ) , differentiating to find the pressure coefficient 
/ Ap = _ 2 80(x,y)' 


ox 


9 * 


and then evaluating the result along the line x = c ^^ . 


Numerical calculations for this base interference effect were carried 
out for the parallel arrangement of bodies to ascertain whether the favor- 
able drag increment would counterbalance the increase in drag previously 
found and indicated in figure J. A rough estimate of the decrement in 
drag was obtained by first plotting the variation in pressure acting on 

the base (illustrated in fig. 8 for a distance parameter 1 = 3.0 and 

several Mach number -sweepback arrangements) obtained by use of the 
formulas in appendix E, and then essentially integrating this pressure 
distribution over the base area affected by the interference flow. The 
decremental drag coefficient as found by this crude approach is 

believed to give a reasonable approximation to the magnitude of the base 

0 £Cr> 

interference effect. Values of the decrement found 

^max^Cmax^ 

by the above-outlined procedure, for the cases presented in figure 8 
were 0.13, 0 .14, and 0.013 for the parameter 0 cot A = 0.10, 0.25, 
and 0.4-0, respectively. Subtraction of these drag decrements from the 
corresponding ordinates of figure 7 result in values that fall beneath 

the "no interference, 4- _» oo" curve presented therein. Additional cal- 
culations covering the range of sweepback, Mach number, and distance- 
be tween-bodies parameter indicated the same result, namely, that the 
base interference effect was generally of sufficient magnitude to over- 
balance the adverse interference effect previously found (see fig. 7) 
and, thus, the overall interference effect on the drag was favorable. Anal- 
ogous calculations for the canted-body arrangement can be carried out with 
the use of the appropriate formulas of appendix E and the procedure previ- 
ously Indicated for the parallel arrangement of bodies. It is apparent that 
the magnitude of the opening angle is a critical parameter with regard to 
the net interference -drag contribution. 


CONCLUDING REMARKS 


An analysis, based on the linearized thin -airfoil theory for super- 
sonic speeds, of the wave drag at zero lift has been carried out for a 
simple two -body arrangement consisting of two wedgelike surfaces, each 
with a rhombic lateral cross section and emanating from a common apex. 


*1 
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Such an arrangement could be used as two stores, either embedded within 
or mounted below a wing, or as auxiliary bodies wherein the upper halves 
could be used as stores and the lower halves for bomb or missile purposes. 
The complete range of supersonic Mach nu mber s has been considered and it 
was found that by canting or orienting the axes of the bodies relative to 
each other a given volume may be redistributed in a manner which enables 
the wave drag to be reduced at the lower supersonic speeds. For purposes 
of comparison, analogous drag calculations for the case of two parallel 
winglike bodies with the same cross-sectional shapes as the canted 
arrangement have been included. Sane consideration has also been given 
to the problem of estimating the favorable (dragwise) interference pres- 
sures acting on the blunt bases of both configurations. In the case of 
the parallel bodies for which calculations were made, this base effect 
seemed more than sufficient to cancel the unfavorable interference on 
the forward part of the configuration. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., July 18, 1957. 
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APPENDIX A 

SUMMARY OF VELOCITY POTENTIALS 
FOR CANTED ARRANGEMENT 

The required velocity potentials for the canted arrangement of 
bodies (see table I) may be conveniently expressed in the following form: 

01A = F I + f 2 + 

0ib = F l + f 3 + F ^ 

01C = £5 + F 2 + 

02A = F 6 + F 2 + F if- 

02B = f 6 + F 3 + 

h = T 1 

K * f 7 + y 8 

where the functions F]_ to Fg are defined as follows: 
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rti/tan* 1 A - (3 


— (x + y tan A)cosh"^ x ,' l:arL — + — y . + 

7 ~ 2 [ P(x + y tan A) 


(x - y tan A) cosh" 1 ?. ten A - ^ 

(3 (x - y tan A) 

. 


rti/cot 2 0 - p 2 


(x + y cot 0)cosh 


(x - y cot eJcosA-l * ,° 0t 9 - 

p(y cot 0 - x) 

, osh -l x cot 9 + ^ 
p(y cot 0 + x) 


= -- j 1 —- — 2 - (kx - y)cos“ 1 ' p2ky r + (kx + yjcos" 1 ~ ff 2ky 
*J&rrL\- P(kX “ y) 0(kx + y] 


Ffi = ^ 2 ) (kx - y)cosh“ 1 ^ + (kx + y)cosh‘ 

jit^l p(y - kx) 


F 7 


(y - x cot A) 


P 2 /cot 2 A - 1 


■ ^ ~ (kx - 


/p 2 k 2 - l 



r 


4 . 


t 


* 
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APPENDIX B 

SUMMARY OF EQUATIONS FOR THE WAVE-DRAG COEFFICIENT 
OF THE CANTED ARRANGEMENT 


For the case wherein all the edges are subsonic, that is, 0 cot A < 1, the formula for the 
wave -drag coefficient is 


pep 


J (p cot A + p tan B) g f" cot h -i £ - (P g cat 2 A+pcot A p tan 8) 

1 , . I n I B qot A « fl tfln 0 


C* /„ \2 XB COt A . . 

V '■max/ '-max; ^(p cot A - p ton 9ty4 - (p tan 9 + p oot A)* L 


p cot A(p tan 9 +■ ?p cot a) cosh _x £ + (p 2 cot 2 A+ P tan 9 p oot A) _ ^ fl ooflll -i 2 - (p 2 ten 2 6 + g tan 9 g cot A) 

P aot A - p ton 0 p ten 9 + 3P cot A P p oot A - p ton 9 


p ton 9(5P ton 9 + P cot A] . 2 + (p 2 tan 2 9 + P too 9 p cot A) _ 2(p cot At p tan 9) 2 . 1 4_ 


OOflh"^ •** W OTP q P WU llj _ fcyp KWh il. T p Wli O/ O00ll“^ ^ + ^ tAP. ^ P COt A) 2 

P cot A - p van 9 p cot A + 5P tan 9 P cot A - p tan 9 4(p ton 0 + p cot A) 


pg cot 2 A 

(p cot A + p tan 0)(3P cot A + p tan 0) caDll -l 2 + (P 2 cot 2 A + p tan a p oot A) 

(p cot A - P tan Q)Jl - P 2 cot 2 A 

fl »4- A A O — — — — Q 4 .„_ n in ..4. a 

p tun ax — p uraii u p bou q T jp but XI 


(p cot A + p ton 9 Jcoott"^* 


2 - (p g oot 2 A+ p tan 6 p cot A) 
p oot A - p tan 9 


* 2p tan 0 cosh -1 


1 - p tan 9 p cot A 
p cot A - p tan 9 


2p tan B(p tan 9 + p cot A) coflh -l It- P tan 9 p cot A 

p aot A - p tan 9 p cot A + p tan 9 0 cot A - p tan 9 


4p a cot 2 A ^jj-l 1 + P g c°t g A 
2p cot A 


p g tan 2 9 


(p cot A ~ p tan 9 - p 2 tan 2 9 


t r rwt a j. a tan nw*- 1 g - (p e tan 2 9-f p cot A P tana) 
,r r ’ p cot A - p tan 9 


H 

VJl 


(equation continued on next pass) 
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(p cot A + p tan 8)(p cot A + 3P tan 8) ^^-1 2 + (p 2 tan 2 B + p cot A p tan 8 ) _ 2 p cot a cost -1 1 ~ P ton 6 P cot A 
3 cot A - p tan 0 p cot A + 3p tan e p cot A - p tan 6 




2P cot A(p cot A + p tan 6 ) COBil -l 1 1 P tan 6 p cot A 4p g tap 2 9 1 + p 2 tan 2 pj 

P cot A - P tan 0 p cot A + P tan 0 p cot A - p tan 6 ** Pp tan B I + 

2jp cot A coe -1 p cot A + 3 tan 8 co6“^" p tan. 8 - (p cot A + P tan ejcos'- 1 - i(p tan 8 + p cot A)|j- 


(Bl) 


As the speed range progresses to the point where the lea ding edge becomes sonic, that is, 
(3 cot A = .1, equation (Bl) reduces to 


( t aax/ c i)a*) 11 (1 - p tan 


B ° D g* a > 2 . L tan e coah-^S tltaih jtan 6(jg_ tanj_ + ^ Z+J_ tan 9 + 

^ e^-OtaneTl^L 1 - p tan 8 ! + 3 p 


p g tan 2 8 


tan 8 


2tl + P ^ tan B) 2 co3h .i 4 
1 - p tan 8 


+ (p tan 6 + l) 2 
4(P tan 6 + 1) 


B 2 ton 2 e 


(1 - p tan 8 )/x - p 2 tan 2 0 


(1 + p tan 8)coeh“ l (2 + p tan 0) + 


(l + P tan 8)(l + 3P tan d) 1 2 + p tan 6 + p 2 tap 2 6 __ 4p 2 tan 2 9 e0 ^_i 1- *■ p 2 tan 2 p 

1 - B ton 9 ' — * * ” " * 


1 + 5P tan 9 1 - p tan I 


23 tan Q 


2 J^p tan 6 coe - *^ p tan P - (1 + p tan 0)coe"4 1(1 + p tan 9 )j|’ 


(B2) 


For the condition in which the leading edge is supersonic but the remaining edges are sub- 

/. « . . \ 

sonic, that is, 3 cot A > 1, p S t the equation for the wave-drag coefficient 


is 


i y 


»-l 


* l 


r * 
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% * 


« * 


* * 




(0 ton 0 + 0 cot A)‘ 


tarn 8 (0 cot A + 0 tan 6 ) a „ . -1 a - 0 Un 6(0 cot A + 0 tan 0 ] 


2 ^ 0 tan 0-0 cot A 


0 cot A - 0 tan 0 


tan 0 (p cot A + 0 tap e)(50 g tan 2 9 + 40 tan 6 0cotA+0 2 cot 2 A) 2+0 tan 0(0 cot A + 0 tan 


(0 cot A - 0 ton 0 ) 2 


0 oot A + J0 ton 9 


2(0 cat A + 0 tan 0)*> en|(h -l 4j (0 cot A + 0 ton 9 ) 2 ^ 0 g cot 2 A 0 2 ton 2 9 - 0 2 cot 2 A 2 - 0 cot a( 0 cot A + 0 tan 8 ) 

{0 cot A - 0 tan 9 ) 2 4(0 cot A + 0 Ian 0 ) J ' Lz cgt 2 A - 1 |(P e °t A - 0 ton 9 ) 2 ” 0 cot A - 0 tan 0 


(P cot A + p tan 9 ) (%? apt A + P tan 9 ) C0B -X 2 + P cot A(p cot A + p t#n 0) + 


(0 cot A - 0 tan 0 ) 2 


tan 8 _i 1-0 cot A 0 ton 6 


50 oot A + 0 ton 0 0 cot A - 0 tan 0 0 cot A - 0 ton 0 


20 tan e(0 cot A + 0 tan 0 } 
(0 cot A - 0 tan 8) 2 


- 

i 1+0 oot A 6 tan 0 

B + — 

P cot A + p tan 0 /“ 


tan2 Q 3 not A + p tan 6 j 2 - p tan e (p cot A P tan 0 ) 
^2 ^^2 0 0 cot A - 0 tan 0 0 cot A - 0 tan 0 


(0 cot A + 0 tan B)(0 cot A + 50 ton 9) h „l 2+0 tan 8(0 cot A+ 0 tan B) 40 2 ton 2 0 


(0 oot A - 0 tan e ) 2 


0 cot A +50 tan 0 (0 cot A - 0 ton 0 ) 2 


_1 1 + 0 2 ton 2 8 
2 p tan B 


2(0 cot A + p tan 0)oos”^ & QQ ^ A * S + gp tan 0 cot p tan 0 


When the ridge line becomes sonic , that is, ^ * co ^ ^ = 1, equation (B3) reduces to 


BCyj 2 I A J" «j P - a cot A 

(tmox/caax ) 2 ’ <t0 cot A |^ a cot 2 A ~3 aot A + 0 cot A - 1 


P - R ftiTh A T ft onrh A - 1_ ft 0nt A _ P i 1 u. Pft r>r\+ A 


+ 1 CO# " 

2 (0 cot A - 1)2 


, o 2 -/i +2 a "1 


( 2 - P ^ A ?^ cosh ' 1 ? - jtLggfr A + P g ” JtS A + 2 ( 2-0 cot A)ca»- 1 (2 - 0 oot A) 

(P aot A - l)^J-3 + 4p cot A - P^ cot 2 A ^ A 
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For the condition in which, both the leading edge and ridge lines are supersonic but the j co 

innermost edge is subsonic, that is, ft ^ * CQ ^ > 1 , p tan 0 < 1, the formula for the ' 

wave -drag coefficient is 
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M 


P 


When the innermost edge 1b sonic, that is, p tan 0=1, equation (B5) reduces to 


P C D 


("^max/ c max) ' 


= 20 cot A 


(p cot A + l) 2 


jp 2 cot 2 A - 1 (p cot A - l)p 2 cot 2 A^^~ cot~A^f _ T) 2_ ^-~^r 


(B6) 


For the speed condition where all edges are supersonic, that is, 0 tan 0 > 1, the drag 
coefficient is given by 


0Cd 


(%a x/ c nmx.y 


2p cot A 


(B tan 9 + B cot A) 2 fl tan 6 


(p 2 cot 2 A)(p cot A - B tan Q)\j(p tan 0 + 0 cot A) 2 - 4 ^0 2 cot 2 A - 1 


(B?) 


It is to be noted that the drag coefficient given by equations (Bl) to (Ff) is based on an area 
which includes the "cutout" portionj that is, the reference area is that of a delta wing having 
the same leading -edge sweepbaek. 


For the special case of 0=0, considerable simplification is introduced into the equations, 
and the formulas for wave-drag coefficient for this case as well as those for the single -vedge 
delta used for comparison purposes in figure 4 are summarized In the following table: 


& 
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APPENDIX C 

SCME SPECIFIC INTEGRALS OF INTEREST OCCURRING 
IN THE DRAG ANALYSIS 


Consider the indefinite integrals 


*1 



-1 cy + d 
ay + b 


<ay 


1 2 - f Josh " 1 SL+J: ay 

U ay + b 

1 3 = fy cos" 1 ■ - - b dy 

U ay + b 

14 = fy cosh -1 cy + d dy 

u ay + b 


where a, b, c, and d are constants subject to the condition that 
ay + b > 0. (Note that this condition is not restrictive inan mnr»h as the 
term can always be made positive by reversing signs in both numerator and 
denominator. ) 

For a 2 > c 2 : 


It = gy - + b cos-1 - cy + - d + - ad - cosh" 1 y(c 2 - a 2 ) + cd - ab * 

a ay + b a i bc - ad l 

1 2 = a y 1 b cosh" 1 cy + d + ad ~ bc cog-1 y(c 2 - a 2 ) + cd - ab 

ay + b a p—z |bc - ad| 

1 3 = g- 2y i- - b2 cos" 1 gy± d - ■ ad - be 

2a2 ay+b 2a(a2 _ c 2) V 

c (a 2 d 2 - b2p2) - 2a2b(a d - be) y(c2 - a 2) + ( c d - ab)* 

2a2( a 2 . c 2)5/2 1*“^ 
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14 


Lj b£ cosh- 1 sr~ + . d 

2a 2 ay + b 


ad - be 
2a(a 2 - c2) 


■JT + 


c(a2d2 _ b2c2) _ 2a^b(ad - be) q y(c 2 
2a 2 (a 2 - c 2 ) 3/2 


a 2) + ( c a. _ a b) 
|bc - ad| 


p p 

For a < c : 


X _ ay+_b oos -l cyj_i + 
ay + b 


ad - be -1 y(c 2 - a 2 ) + (cd - ab) 

cos — r— - — — — - 

~2 |bc - ad| 


-1 y(c 2 - a 2 ) + (cd - ab)* 


_ ay + b , cy + d ad - be , . 

Ip = cosh — + — - - cosh" , . 

a ay + b a ^ c 2 _ a 2 l bc - ad | 


1, = a2 y 2 -.l 2 . cos- 1 + - ad -_ bc .. J~Y 

' 2a 2 ay + b 2a(c 2 - a 2 ) ’ 


c(a 2 d 2 - b 2 c 2 ) - 2a2b(ad - be) COB -l y(c 2 - a 2 ) + cd - ab 


2a 2 (c 2 - a 2 ) 


3/2 


be - ad| 


14 


= s2y2 - b2 cosh- 1 ° y -t a + - £ d - b ? -R . 
2a 2 a y + b 2a(c 2 - a 2 ) 


c(a 2 d 2 - b 2 c 2 ) - 2a 2 b(ad - be) ^g^-l y(c 2 - a 2 ) + cd - ab^ 
2a 2 (c 2 _ a 2)3/ 2 C0S l lc " ad l _ 


where the quantity Y = (c 2 - a 2 )y 2 + 2(cd - ab)y + d 2 - b 2 and the 
asterisk indicates that if the inverse hyperbolic coBine should be 
required of a negative number -N then the following relationship must 
be used: 


cosh -1 (-N) = -cosh -1 N 
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APPENDIX D 

SUMMARY OF FORMULAS FOR VELOCITY POTENTIALS AND WAVE 
DRAG FOR THE PARALLEL ARRANGEMENT 


The basic potentials required to evaluate the wave drag of the 
parallel arrangement of surfaces (see fig. 3) nay be obtained for the 
case of subsonic leading edges involving interference effects by the 
procedure outlined in the body of the paper. At all supersonic Mach 
numbers for which the particular region of one surface under considera- 
tion is outside the domain of interference from the opposite surface , 
the potential expressions may be readily reduced from reference 3 . For 
convenience, all the applicable potential expressions are summarized in 
the following formulas: 

Subsonic leading edge; region outside interference field: 


0 _ . ^(^max/ c max) 


2^1 - p 


cot^ A L 


(y - x cot Alcosh-l V -ye 2 cot A , - 

P(x cot A - y) 


(y + x cot A)cosh-l V ^ cot A 


p(x cot A + y) 


(Dl) 


Subsonic leading edge; region within interference field: 


0 _ X. Ofrjax/ c max) 


2rt^L - p2 cot 2 A 


(y - x cot AjcostT 1 gotA ' - 

P(x cot A - y) 


(y + x cot A) cosh" ^ x , + — co ' t: ' A + 

p(x cot A + y) 

(y + l - x cot AjcosH- 1 x (y t DP 2 C °* A . 

p(y + l - x cot A) 


(y + l + x cot A)cosh“- L x+ ^ + cot A 


p(y + l + x cot A) 


(D2) 
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Supersonic leading edge; region outside Mach traces from apex: 


0 = v (w/w) fr - x cot A ) 

2^2 cot2 A - 1 


(D3) 


Supersonic leading edge; region within Mach traces from apex: 


0 _ ^ C^max/ Siax ) 

2n^p 2 cot^ A - 1 


(y - x cot A)cos -1 * ~ ° ot A 

3(x cot A - y} 


f . . \ _i x + yS^ cot A 

(y + x cot A) cos x r 

P(x cot A + y) 


(D4) 


The formulas for wave -drag coefficient are based oh the plan-form 
area of both surfaces and are given as follows: 


For 0 < p cot A 5 


< 1 . 


1+1 


7 - C — ... .. = . 1 . B cyt A | . 

(W/craax) 2 * y/l - p2 oot 2 A L 


2 cosh 


- 1.1 + 


*• P 2 cot£ A _ (-L + if coBh- 1 1 + (a + X ) p2 COtg A 
2g cot A \2s ) (I+2)peofA 


4s 


4(1 - p 2 cot 2 A) + i(l + P 2 cot 2 A) i 1 + (i - l)p 2 cot 2 A 

B cosh — - 


1 - P 2 cot 2 A 


— B cot A 
a 


4(1 - p 2 cot 2 A) - 1(1 + p 2 cot 2 A) 1 - (1 + l)p 2 cot 2 - A 

r* - cosh'l ^ >- 

l _ b 2 2 » 


1 p cot A 


f 1 \2 -1 1 ' (s ‘ 00t2 A 

f_L _ i) cosh 1 -15 1 + 

\2s / (1 - 2^p cot A 


cot A 


h , # 3 ( f ^ W ^- F 7 ¥ W z **) - 

yl - p 2 cot 2 A'' J 

2 sin -1 p cot A + cos -1 j^i + if cot aJ- cos -1 j^l - if cot i\jl 


(D5) 
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For 



< (3 cot A 



m _ 1 

( W r- j C ttw v ^ * 


9 cot A 

I 4(1 - g 2 cot 2 A) + g^l + g 2 cot 2 A) 

^1 - g 2 cot 2 A 

1+6 1 - g 2 cot 2 A 


cosh" 


^ p cot . 
s K 


!l - 1 - (& - A _ 

P 1 + e 2 “t 2 A 

'* 7 (M 

g cot A 

2g cot A 

i i 1 ~ (h - l ) 2 ? 2 °° t£ A 

- 2 sin -1 g 

cot A - cos-|i - ljg ( 

- g 2 cot 2 A 


For 


1 - 1 


< p cot A < 1: 


(3C D 


= 2, 


('tnax/'Tnax) 2 * 


■fi-SS* - 4 . . cosh- 1 fi i cot . 2 A + 1 + sin- 1 p 

2p cot A 


/l - p 2 cot 2 A 


For p cot A ^ 1: 


PC D 


(^mx/ c max) i 


= 1 


cot 2 A 

i 

(d6) 

' 

cot A? 

«* 

(D7) 

(D8) 
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APPENDIX E 

EXPRESSIONS FOR INTERFERENCE 
PRESSURES ACTING ON RASE 


The formulas in this appendix give the interference pressures (in 
coefficient form) acting along the base line in the plane of synmetry, 
that is, x = c ^ y and z = 0. 


Canted Arrangement 

For the .case where all the edges are subsonic, that is, p cot A < 1, 
the formula is 


(tmaac/cmaac)^ 11 Xaa© int 


Ji 

(P cot A + p tan 0) 2 

x(ji cot A - p tan 8) 

^4 - (p cot A + p tan S ) s 


3 2 cot 2 A 




„-l 1 + (a)S 2 “t 2 A B 2 tan 2 0 . _i 

- coah A — - - — ° — coah 1 - , 

p 2 cot 2 A (l + jQp tot A ^L-p 2 ” tai^0 P tan 0 + 


2 + (p tan 6 + p cot A)(j[)p cot A 
p tan 8 + (l + 2 £)p cot A 

1 + (f)(0 tot A)p tan 0 


cot A 


(SI) 


When the leading edge becomes sonic, that is, p cot A = 1, equa- 
tion (El) reduces to 


KfcA 


2 

[ q t 8 tm 8) 2 


coah-1 2 +»»“»> 

Int «(1 - 0 tan 8 ) 

^\j( 1 - 3 tern 8) (3 3 tan 8 ) 

l+ptan 8 + 2 ^ 

p 2 tan 2 8 v-1 1 + jc 0 tan 8 

1 -S 


^1 - p 2 tan 2 8 

0 tan 0 + £ J 

i+ fJ 



(E2) 


For the condition where the leading edge is supersonic but the 
remaining edges are subsonic, that is, p cot A > 1, p ( ^ n . . Q ± PPji , < 1 ; 

the interference pressure is given by 


6 


2 

(0 cot A + 0 tan 8) 2 „. k -l 

2 + (f)(P cot A)(P tan 8 + 0 cot A) 

(tjjiax/cmax) \ ^ /haso int 

«(3 cot A - 3 tan 0) 

^4 - (p cot A + p tan e)" 2 

0 tan 8 + (l + S^0 cot A 



P 2 cot 2 A 1-+ (b)^ 2 cot£ A B 2 tan 2 8 

_i 1 + (|)(P cot A)0 tan 8 



yp 2 cot 2 A - 1 ! 

fl + |^P cot A _ p 2 tan 2 8 

3 tan 8 + cot A 


(E5) 



NACA TN 4-120 


27 


When the ridge line becomes sonic, that is, 
equation (E3) reduces to 


r tan 9 + cot A 

v 5 


)■ 


CW?w)(t) 


2 

Li 

n 

;Dp cot A b 2 cot 2 A „ 

0,-1 1 + ($>* C ° t£ A 

base lnt *(P cot A - p tan 8) 

J 

i +< 

'|)pcotA ^2 CQt g A - 1 

(l + j^p cot A 


t 2 ten 2 a co6h -l cot A)p ten 8 

£ - p 2 ten 2 6 P tan 0 + (£)p cot A 


0*) 


For the condition in which both the leading edge and ridge lines 
are supersonic but the innermost edge is subsonic, that is, 
p ^tan 8 + c _ ot > p and J3 tan 8 < 1, the pressure formula Is 


p 

( [%auLj c tmx ) 


\l] 


base int 


2 

(p cot A + P tan 0)^ 

«(p cot A - p tan ej 

^(P cot A + p tan 0) 2 - k 


2 + 


(a'X^ CQt tan ft + p cot A) 

p tan 6 + ^1 + 2 gjp cot A 


pg cot 2 A CM -1 1 ♦ (£)p 2 cot g A _ pg tan 2 8 coall .i 1 + (£)(P CQt A)p tan 6 
y oot 2 A - 1 CM (l + |)P cot A h . pg tan 2 8 P tan 8 + (^)p cot A 


At the higher Mach numbers (where the inner edge becomes sonic or 
supersonic) the interference field does not affect the base region of 
the opposite panel. 


Parallel Arrangement 

For the parallel arrangement of bodies the expression for the 


interference 

pressure is 



e m 

p cot A 

. 1 - (1 + j)p 2 cot 2 A 

OOS.-1 

( t totx/ c mx)V‘l/baee Int _ p 2 ^2 A 

(Z. t l-l)pcotA 

(1 + i * l )p cot A 


(B6) 
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TABLE I.- SUMMARY OF CASES TREATED, MATHEMATICAL CONDITIONS, 
AND APPLICABLE POTENTIAL EXPRESSIONS 


B 

e 

Description 

Mathematical conditions 

Applicable 

potential 

expressions 

i 

All edges subsonic 

B cot A < 1 

01B> 02B 

ii 

Leading edge supersonic 
but other edges subsonic 

P tan 0 + p cot A < 2 
P cot A > 1 

^lA* ^2A' ^3 

hi 

Only inner edge subsonic 

P tan 9 < 1 
P tan 0 + p cot A > 2 

01C> ^3> 

IV 

All edges supersonic 

P tan 0 > 1 
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Figure 3 . - Pertinent geometry for parallel arrangement of surfaces and sketches indicating 

various combinations possible for arbitrary values of Mach number, aspect ratio, and distance 
parameter Z/s. 
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Figure 5,- Variation of wave -drag coefficient with Mach number -sweepback 
parameter for canted arrangement of surfaces. 
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Figure 8.- Illustrative variations of the interference pressures acting along the base of 
one panel due to the presence of the opposite panel for the parallel arrangement of 
bodies, i/s = 3-0; z ° 0; x = c max* 






